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'  The  MILLS  105  nr.  main  tank  gun  transmitter  was  tested  for  the  potential 
of  using  the  system  as  a  marksmanship  trainer.  The  tests  were  carriec.  out 
over  a  3000  meter  range  at  the  Space  Shuttle  runway  at  Kennedy  Space  Center. 
Experimental  evidence  collected  in  CV  laser  experiments  snow  that  the 
statistical  distribution  used  in  the  MILLS  development ,  i.e.,  Log-Normal, 
is  valid  only  over  short  ranges.  At  long  ranges  the  Negative  Exponential 
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20.  Abstract  (cont’d) 


.  distribution  was  measured.  A  new  universal  matn  model  was  developed 
which  accurately  aescribed  the  statistics  of  the  fluctuation  laser  signal 
along  the  entire  range  up  to  3000  meters.  Using  this  new  universal  model 
the  cunmulative  fading  was  calculated.  This  analysis  showed  that  it  is 
the  intermediate  ranges  which  have  the  greatest  fading.  Based  upon  the 
experimental  eviaence  and  the  math  model  we  conclude  that  the  MILLS  sysie 
is  not  sufficiently  accurate  for  long  range  gunnery  marksmanship  training, 
but  that  a  new  C.V  laser  transmitter  and  receiver,  developed  by  kTLC  in  a 
parallel  project,  should  meet  the  required  accuracy. 
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SUMMARY 


This  document  reports  on  a  research  project  directed  by  PM  TRADi: 
and  performed  jointly  by  the  Engineering  and  Industrial  Experiment 
Station  of  the  University  of  Central  Florida  (UCF)  and  the  Research 
Department  of  the  Naval  Training  Equipment  Center  (NTEC).  The 
project  centered  on  the  study  of  the  use  of  the  MILES  system  as  a 
marksmanship  trainer  for  long-range  tank  and  anti-tank  gunnery. 

The  study  evaluated  the  MILES  105mm  laser  transmitter  with  the 
MILES  105mm  standard  target.  During  the  tests  the  atmospheric  optical 
channel  was  simultaneously  monitored  by  a  CW  laser  propagating  parallel 
to  MILES  105mm  beam.  The  test  results  show  the  MILES  beam  is  much 
wider  than  predicted  by  the  analysis  in  the  document  MILES  A002,  June 
1980.  Also  a  comparison  of  the  UCF  and  the  NTEC  test  results  show 
that  the  background  sun  illumination  significantly  reduces  the  effective 
range. 

The  study  also  evaluated  the  statistical  model  (Log-Normal  distribution) 
used  as  a  basis  of  the  MILES  design  and  for  prediction  of  its  operational 
characteristic.  The  evidence  collected  during  experiments  conducted  at 
the  Space  Shuttle  runway  at  Kennedy  Space  Center  shows  that  the 
Log-Normal  distribution  was  valid  only  up  to  ranges  of  about  180  meters. 

At  long  ranges,  under  conditions  of  strong  atmospheric  turbulence,  the 
Negative  Exponential  distribution  was  measured.  The  UCF  team  developed 
a  universal  math  model  which  agreed  with  both  the  distribution  at  the 
short  and  long  ranges  as  well  as  the  measured  distribution  at  inter¬ 
mediate  ranges.  The  cumulative  fade  was  then  calculated  from  the 
universal  math  model.  The  results  show  that  the  intermediate  ranges 
have  the  largest  fading. 
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The  recommendations  for  MILES  add-ons  are:  to  replace  the  MILES 
detectors  with  photo-diodes,  to  replace  the  transmitter  with  a  high 
pulse  rate  CW  laser  diode  transmitter,  and  to  replace  the  receiver  with 
a  corresponding  phase  locked  loop  or  matched  filter  receiver.  These 
systems  could  be  made  completely  compatible  with  the  MILES  encoder 
and  decoder.  The  threshold  of  this  receiver  and  the  pulse  rate  can  be 
set  by  analysis  of  the  UCF  universal  math  model  for  the  laser  propaga¬ 
tion  communication  channel.  The  study's  results  support  the  conclusion 
that  the  current  MILES  system  is  not  sufficiently  accurate  to  be  used  as 
a  marksmanship  gunnery  trainer.  The  analysis  supports  the  contention 
that  the  required  accuracy  can  be  obtained  by  using  a  high  pulse  rate 
CW  laser  diode  transmitter  and  the  appropriate  CW  receiver  as  proposed 
by  the  NTEC  research  team. 

The  PM  TRADE  project  director  was  A.  J.  Boudreaux.  The  NTEC 
research  team  leader  was  A.  H.  Marshall,  and  the  UCF  team  head  was 
R.  L.  Phillips. 
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SECTION  I 


INTRODUCTION 


This  document  reports  on  a  project  which  undertook  the  evaluation 
of  the  MILES  system  as  a  long-range  marksmanship  gunnery  trainer. 

The  project  consisted  of  measurements  on  a  MILES  tank  main  gun  trans¬ 
mitter,  analysis  of  the  mathematical  models  of  the  laser  propagation, 
analysis  of  the  MILES  receiver  binary  union  decoder,  and  the  develop¬ 
ment  of  a  new  CW  laser  transmitter  and  receiver  which  could  be  used  as 
an  add-on  for  marksmanship  uses  of  MILES.  The  project  was  organized 
in  two  research  teams.  The  University  of  Central  Florida  (UCF)  team 
performed  the  mathematical  analysis  and  the  experimental  verification  of 
the  mathematical  models  of  the  laser  propagation  phenomenon  upon  which 
much  of  the  MILES  system  design  is  based.  The  Research  Department 
of  the  Naval  Training  Equipment  Center  (NTEC)  conducted  the  develop¬ 
ment  of  the  CW  laser  transmitter  and  receiver  hardware.  In  addition 
both  teams  took  extensive  measurements  on  the  MILES  tank  main  gun 
simulator. 

The  MILES  system  was  designed  and  constructed  as  a  short  range 
tactical  laser  weapon  fire  simulator.  It  was  designed  not  to  simulate  the 
flight  of  a  particular  round,  but  rather,  to  simulate  the  probability  of  a 
kill  in  an  engagement.  The  "MILES  kill"  may  occur  because  of  a  success¬ 
ful  detection  and  decode  of  a  MILES  code  word  by  any  of  the  MILES 
detectors  on  the  target.  To  ensure  kills  in  a  tactical  engagement,  the 
transmitted  laser  beam  was  made  sufficiently  wide  to  make  the  kill  zone 
actually  wider  than  the  target.  In  other  words,  the  weapon  can  be 
aimed  at  a  point  actually  off  the  target  and  a  kill  could  still  be  obtained. 
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If  MILES  is  to  be  used  in  gunnery  marksmanship  training,  a  much 
smaller  and  well-defined  kill  zone  will  be  required.  The  kill  zone  is 
defined  by  an  inter-play  between  the  transmitter  and  receiver  character¬ 
istics,  but  it  is  the  receiver  structure  which  has  the  dominant  effect. 

The  basic  receiver  can  be  divided  into  three  portions  —  the  detector,  the 
receiver,  and  the  decoder  which  is  used  to  decipher  the  MILES  codes. 
The  mathematical  model  used  to  describe  the  laser  atmospheric  channel 
which  served  as  a  basis  for  their  design  and  performance  evaluation  of 
the  decoder  and  the  codes,  was  the  Log-Normal  statistic.  If  any  MILLS 
add-ons  were  to  be  designed,  the  validity  of  this  model  would  require 
testing.  Some  of  the  experiments  carried  out  to  investigate  the  statis¬ 
tics  of  long-range  laser  propagation  were  described  in  a  previous  report 
[1].  Those  experiments  were  augmented  by  further  experiments  con¬ 
ducted  in  August  1980  at  Kennedy  Space  Center  Shuttle  Landing  Facility 
The  results  of  those  investigations  are  reported  in  Section  III  of  this 
document. 

The  essence  of  the  long-range  detection  and  decoding  problem  is 
the  determination  of  transmitted  MILES  code  from  a  random-fading 
received  signal  in  the  midst  of  receiver  noise.  Two  techniques  widely 
used  in  communication  systems  to  extract  a  fading  signal  in  the  pres¬ 
ence  of  noise  are  phase-locked  loops  and  matched  active  filters.  These 
two  approaches  were  investigated  by  the  research  team  at  NTEC.  Both 
of  these  receiver  types  require  a  CW  laser  transmitter  operating  at  a 
high  pulse  rate.  The  required  detector  for  both  of  these  techniques  is 
a  high  speed  device  such  as  a  reversed  biased  photo-diode  or  photo¬ 
avalanche  diode.  The  parameters  that  must  be  determined,  for  either 
the  phase-locked  loop  or  the  matched-active  filter,  are  the  receiver 
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threshold  and  the  transmitter  pulse  repetition  frequency  (PRF).  If 
these  systems  are  to  be  compatible  with  the  MILES  encoder  and  decoder, 
the  systems  must  operate  at  a  multiple  of  the  MILES  PRF  of  3  KHz. 

For  instance,  if  twenty  pulses  were  required  by  the  receiver  to  acquire 
the  signal,  the  necessary  PRF  would  be  60  KHz.  Then  every  time  the 
MILES  encoder  commanded  the  transmitter  to  send  a  pulse,  correspond¬ 
ing  to  a  MILES  code  "one",  the  CW  transmitter  would  transmit  twenty 
pulses  during  the  same  period.  The  detector  would  then  detect  the 
signal  and  the  receiver  would  lock  onto  the  twenty  pulses  and  then 
pass  a  single  pulse  onto  the  MILES  decoder  at  the  MILES  PRF  of  3  KHz. 
Hence  the  new  CW  laser  transmitter  and  detector/receiver  would  be 
"invisible"  to  the  MILES  encoder  and  decoder. 

So  as  to  determine  the  best  PRF  for  the  transmitter,  the  optimum 
receiver  threshold  and  the  required  transmitter  beam  width,  a  mathemat 
ical  model  is  required  to  accurately  predict  the  statistical  scattering  of 
the  atmosphere.  The  required  models  are  the  probability  density  func¬ 
tion,  the  cummulative  fading,  and  the  average  fade  period  length.  The 
models  must  be  universally  accurate  at  all  extended  ranges  if  the  system 
is  to  simulate  the  flight  of  a  particular  round  as  required  in  marksmansh 
gunnery  training. 

Section  II  reports  on  the  measurements  made  on  the  MILES  main 
tank  gun  simulator.  The  atmospheric  turbulence  conditions  and  the 
background  sun  were  simultaneously  monitored  and  their  effects  on  the 
MILES  simulator  are  discussed. 

Section  III  develops  the  universal  mathematical  model  for  laser 
beam  propagation  over  ranges  up  to  3  Km.  This  section  also  discusses 
the  fading  time  of  the  scattered  laser  signal. 
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SECTION  II 


MEASUREMENTS  OF  THE  MILES  SYSTEM 

Measurements  were  made  on  a  MILES  105  mm  tank  main  gun  laser 
transmitter  and  an  armored  personnel  carrier  (APC)  receiver  system. 
Data  was  taken  on  the  system  on  a  darkened  indoor  range  and  on  a 
long  outdoor  range. 


EXPERIMENTAL  LAYOUT 

The  indoor  tests  were  made  at  a  range  of  exactly  100  feet  in  com¬ 
plete  darkness.  The  MILES  105  mm  transmitter  was  rigidly  mounted  on  a 
tripod  one  meter  above  the  floor.  The  target  was  the  MILES  APC  receiver 
using  only  one  detector.  The  other  detectors  were  all  taped  to  prevent 
any  light  leak  whatsoever.  Initially,  the  transmitter  was  aimed  directly 
at  the  detector  and  then  moved  off  center  both  horizontally  and  verti¬ 
cally  in  one-inch  steps.  At  each  position  the  transmitter  was  fired  ten 
times.  If  all  ten  shots  were  kills,  the  transmitter  aim  point  was  moved 
another  inch  away  from  the  detector.  When  one  miss  was  recorded,  the 
location  was  noted.  The  transmitter  aim  point  was  moved  in  each  of  the 
four  directions  off  of  the  detector  as  shown  in  Figure  2-1.  These 
measurements  then  established  the  divergence  of  the  90%  kill  zone  in  the 
near  field  of  the  transmitter. 

The  long-range  tests  were  carried  out  at  the  Space  Shuttle  runway 
at  NASA's  Kennedy  Space  Center.  The  runway  is  16,000  feet  long  and 
150  feet  wide  and  completely  level.  It  is  constructed  of  concrete  with 
grooves  cut  across  it  every  inch  for  the  entire  length.  Hence,  the 
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exact  range  over  which  the  tests  were  made  was  measured  very  accu¬ 
rately.  The  area  around  the  runway  was  cleared  of  all  structures  and 
all  vegetation,  except  grass,  for  distances  of  several  hundred  feet. 

These  features  of  the  runway  allowed  the  MILES  tests  to  be  conducted 
over  a  completely  homogeneous  range  ensuring  the  test  results  would  be 
repeatable.  A  schematic  of  the  Space  Shuttle  runway  is  shown  in 
Figure  2-2. 

The  MILES  105  mm  laser  transmitter  was  mounted  on  an  optical  rail 
along  side  of  a  HeNe  (wavelength  0.6328  mm)  CW  laser  also  mounted  on 
an  optical  rail.  The  two  rails  were  then  attached  to  the  head  of  a  tripod 
which  could  be  rotated  both  in  the  horizontal  and  vertical  planes  with 
an  accuracy  of  an  arc  minute.  The  lasers  were  pointed  down  the 
length  of  the  runway  at  a  height  of  1.3  meters  above  and  level  with  the 
surface. 

The  MILES  APC  receiver  belt  system  was  mounted  on  a  standard 
MILES  105  mm  target  as  described  in  the  Xerox  CDRL  Item  A002,  Trainer 
Engineering  Report  for  MILES  Volume  I,  June  4,  1980.  The  APC  receiver 
belt  was  attached  to  the  target  with  five  of  the  six  detectors  facing  up 
range,  and  the  sixth  detector  completely  taped  as  prescribed  in  Xerox 
A002.  The  entire  target  was  then  mounted  vertically  on  the  side  of  a 
van  at  a  height  of  three  feet  above  the  runway  surface.  The  target 
mounted  on  the  van  is  shown  in  Figure  2-3.  The  target  was  moved  in 
one-half  meter  increments  across  the  MILES  105  mm  transmitter  beam 
and  at  each  new  position  ten  shots  were  fired.  This  procedure  was 
repeated  until  one  shot  was  missed  out  of  ten  shots  fired.  These 
measurements  were  performed  across  the  beam  thereby  defining  the  90% 
kill  zone  of  the  beam  at  that  range.  The  measurements  were  taken  at 
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500  meter  increments  starting  at  500  meters  and  continuing  down  range 
to  3000  meters. 

During  the  MILES  measurements,  the  background  solar  illumination 
was  monitored.  Also  monitored  were  the  temperature,  wind  speed,  and 
wind  direction. 

The  CW  laser  system  parallel  to  the  MILES  beam  allowed  for  simul¬ 
taneous  measurements  of  the  atmospheric  channel  through  which  the 
MILES  signal  was  being  transmitted.  The  CW  laser  signal  was  measured 
at  a  short-range  station  and  at  a  long-range  station.  The  short-range 
station  was  located  at  a  fixed  position  183  meters  from  the  transmitter. 

At  this  station  special  electronics  sampled  the  CW  laser  beam  every  minute 
and  calculated  the  average  signal  level  and  the  variance  of  the  CW  laser 
signal  fluctuations  induced  by  atmospheric  turbulence.  The  long-range 
station  was  always  located  down  range  at  the  same  position  as  the 
MILES  receiver  target.  The  laser  signal  was  monitored  at  this  station 
by  a  MINC-11  computer  built  by  Digital  Equipment  Corporation.  The 
MINC-11  computer  calculated  the  first  five  statistical  moments  of  the 
fluctuating  CW  laser  signal.  Before  the  tests  of  the  MILES  system 
began,  the  long-range  station  was  placed  alongside  of  the  short-range 
station.  Using  the  first  five  statistical  moments,  it  was  established  that 
at  183  meters  from  the  transmitter  the  short-range  station  would  be 
located  within  the  range  for  which  the  Log-Normal  mathematical  model 
for  the  laser  scattering  would  be  valid.  In  this  regime  only  the  mean 
value  and  the  variance  of  the  fluctuating  signal  are  needed  to  calculate 
the  so-called  atmospheric  turbulence  structure  constant,  C*.  This  con¬ 
stant  is  the  measure  of  the  strength  of  the  atmospheric  turbulence, 
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SHORT  RANGE  STATION 


1 


where: 


turbulence 

c2  „-2/3 

n 

weak 

< 

10"15 

moderate 

10‘14 

strong 

> 

io‘12 

Hence  using  the  data  from  the  short-range  station,  the  strength  of  the 

turbulence  present  in  the  optical  channel  can  be  calculated.  The  homo- 

2 

geneity  of  the  Shuttle  runway  assured  the  measured  Cn  at  the  short- 
range  station  was  valid  for  the  entire  range.  The  statistical  moments 
from  the  long-range  station  along  with  the  data  from  the  short-range 
station  allowed  for  the  complete  characterization  of  the  scattering  statis¬ 
tics  of  the  CW  laser  channel  and  thereby  the  statistics  of  the  parallel 
MILES  laser  communication  channel.  The  experimental  layout  of  the 
tests  at  the  Space  Shuttle  runway  are  shown  in  Figure  2-4. 

RESULTS 

The  90%  kill  zone  for  the  short-range  indoor  tests  is  shown  in 
Figure  2-1.  The  near  field  90%  kill  zone  has  an  effective  horizontal 
divergence  of  20  milliradians  and  a  vertical  divergence  of  14  milliradians 
at  100  feet. 

The  long-range  outdoor  measurements  of  the  90%  kill  zone  is  shown 
in  Figure  2-5.  The  dotted  curve  shows  the  predicted  90%  kill  zone  as 
depicted  in  Figure  F-16,  page  F-34,  of  Xerox  A002,  Volume  I,  June  4, 
1980.  The  Xerox  curve  used  for  comparison  was  the  50  Km  visibility 
curve  since  there  were  unlimited  visibility  conditions  during  the  MILES 
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Fig.  ?.-S.  The  soliH  curve  shows  the  907.  prohnhilitv  curve  for  the  I  OS  inm  transmitter  .is  measured 
at  the  Shuttle  runwnv.  Measurements  were  made  on  I  v  out  to  1000  meters  hut  the  data  shows  t  >  i 
the  00?'.  kill  zone  actually  extends  s'ightly  further.  The  dotted  curve  is  the  theoretically 
predicted  curve  from  the  Xerox  A002  document. 


tests  at  the  Shuttle  runway, 
varied  only  slightly  from  10 


2 

The  turbulence  structure  constant  Cn 

-2/3  .13  -2/3 

m  '  to  10  m  which  means  the 


2 

turbulence  remained  very  strong.  The  Cn  was  calculated  using  the 


equation 


2 


Ln  (<I2>/<I>2) 

- Tit - ITTB 

0.5  K  /  L  1 


(2.1) 


where  <I2>/<I>2  is  the  normalized  second  moment  obtained  from  the 
measurements  of  <I2>,  and  <I>  at  the  short  range  station  for  the  CW 
laser  beam,  K  =  2n/A  (A  is  the  wavelength),  and  L  is  the  range  (183 
meters). 

During  the  tests  on  the  MILES  system  the  background  sunlight  was 
monitored.  The  background  measurements  were  made  with  a  40X  OPTO¬ 
METER  from  United  Detector  Technology  Inc.  The  device  was  calibrated 
in  lumens/square  foot  or  footcandles.  Typically,  the  sun  illuminance  on 
the  MILES  target  varied  from  103  to  7xl03  lumens/square  foot.  The 
variation  was  due  to  the  cloud  movements  in  the  partially  overcast  sky. 
The  partially  overcast  sky  and  the  fact  that  the  target  was  facing 
almost  south  (150°  from  north)  meant  that  very  little  of  the  direct 
sunlight  was  incident  on  the  target.  This  should  be  contrasted  with 
the  experimental  layout  of  the  MILES  tests  conducted  by  NTEC.  In 
those  tests  the  MILES  receiver  belt  faced  almost  due  west  [2].  Since 
those  long-range  tests  were  conducted  later  in  the  day,  the  full  sun 
shown  directly  onto  the  target. 

The  effect  of  background  sunlight  on  the  MILES  unbiased  solar  cell 
can  be  severe.  The  effects  of  background  illumination  has  been  studied 
in  an  earlier  publication  of  Mallette  and  Phillips.  A  reprint  of  that 
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publication  is  enclosed  in  Appendix  D.  The  effects  of  temperature  on 
the  MILES  detector  has  also  been  explored  and  reported  in  Reference  16) 
The  effect  of  high  temperature  and  especially  background  illumination 
is  to  reduce  the  sensitivity  of  the  detector.  It  is  this  reduced  sensitivit 
due  to  the  greater  sun  background  light,  which  caused  the  NTEC  test? 
to  show  an  effective  90%  kill  range  of  only  2,500  meters  while  the  UCF 
tests  at  Kennedy  Space  Center  show  an  effective  90%  range  of  over 
3,000  meters  as  shown  in  Figure  2-5.  The  NTEC  experimental  results 
are  shown  in  Figure  2-6,  for  comparison. 

The  UCF  measured  90%  beam  width  shows  several  ’’bulges"  in  the 
beam.  These  "bulges"  seem  to  be  due  to  the  tri-gaussian  nature  of  the 
transmitter  beam  as  described  in  Xerox  A002.  The  evidence  to  support 
this  is  that  the  first  "bulge"  correspnds  to  the  90%  kill  zone  divergence 
as  measured  on  the  100-foot  range.  The  second  bulge  occurring  at 
about  1,500  meters  is  not  as  pronounced  as  the  first  due  to  the  natural 
beam  spreading.  The  "bulge"  at  the  longest  range  corresponds  to  the 
third  gaussian  component  of  the  tri-gaussian  beam. 

DISCUSSION 

The  data  taken  in  the  tests  at  Kennedy  Space  Center  showed  a 
much  greater  range  and  wider  kill  zone  for  the  MILES  system  than  did 
the  NTEC  test  results.  This  is  due  to  the  much  lower  background  sun 
illumination  of  the  MILES  detectors  in  the  UCF  tests. 

The  UCF  tests  showed  that  the  beam  width  varies  along  the  range. 
This  was  due  to  the  inherent  shape  of  the  laser  diode  output  beam  in 
the  MILES  transmitter.  The  beam  width  shows  three  bulges  character¬ 
istic  of  a  tri-gaussian  beam  shape. 


Fig.  2-7.  (a).  CW  laser  transmitter  used  to  monitor  atmospheric  optical 

channel  during  MILES  test.  (b) .  CW  laser  as  seen  from  down  range 
from  MILES  target  location. 
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SECTION  III 


UNIVERSAL  MODEL  FOR  TURBULENCE  LASER  SCATTERING 

An  accurate  universal  (valid  at  all  ranges)  mathematical  model  is  a 
stronger  requirement  for  marksmanship  training  than  for  tactical  train¬ 
ing  due  to  the  small  well  defined  kill  zone  needed  for  the  marksmanship 
trainer.  The  tactical  training  weapon  fire  simulator  requires  a  kill  zone 
that  encompasses  the  entire  target  while  the  marksmanship  trainer 
requires  a  kill  zone  which  is  only  a  small  portion  of  the  target.  Hence, 
if  a  matrix  of  independent  detectors  is  used  on  the  target,  as  sug¬ 
gested  by  Reference  [2],  this  means  that  only  one  detector  would  be 
illuminated  by  any  one  laser  shot.  This  is  contrasted  with  the  MILES 
tactical  system  which  may  have  several  detectors  illuminated  at  one  shot 
and  their  individual  signal  outputs  summed  before  the  threshold  is 
encountered.  The  summed  outputs  of  the  multiple  MILES  detector 
reduces  the  variances  of  the  turbulence-induced  laser  signal.  This 
would  not  be  the  case  for  the  gunnery  marksmanship  target.  Since 
each  detector  would  be  independent  and  have  an  independent  threshold, 
the  variance  of  the  detected  signal  will  be  the  same  as  that  of  the  laser 
light  itself  propagating  through  the  atmospheric  turbulence.  This 
variance  will  be  higher  than  that  in  the  MILES  system. 

Also,  as  pointed  out  in  Section  II  and  in  Reference  [2],  the  un-  j 

biased  detectors  of  the  MILES  system  are  not  adequate  for  the  marks¬ 
manship  training  due  to  their  loss  of  sensitivity  in  the  presence  of 
background  sun  illumination.  To  obviate  this  problem,  reversed  biased 
photo-diodes  were  suggested.  These  small  photo-diodes  are  very  much 
smaller  than  the  MILES  one  square  centimeter  solar  cell  detectors.  As  a  i 
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consequence,  the  averaging  effects  realized  by  the  large  MILES  detec¬ 
tors  which  reduced  the  fading  effects  of  the  turbulence  would  not  occur 
in  the  smaller  photo-diodes.  This  averaging  phenomenon  had  an  addi¬ 
tional  mitigating  effect  on  the  variance  of  the  turbulence-induced  fluctu¬ 
ation  of  the  laser  signal.  It  is  because  of  lack  of  these  averaging 
effects  in  the  required  detector  for  marksmanship  that  the  electrical 
signal  after  detection  will  be  identical  to  the  statistics  of  the  turbulence- 
induced  laser  light  fluctuations. 

In  a  previous  report  [1]  a  universal  mathematical  model  was  con¬ 
structed  for  laser  propagation  through  turbulence  which  could  be  used 
for  the  statistical  analysis  of  laser  weapon  fire  simulators.  That  mode! 
was  based  upon  data  gathered  during  March  1980  through  experiments 
conducted  at  Kennedy  Space  Center's  Shuttle  Landing  Facility.  Although 
the  universal  model  developed  in  the  earlier  report  did  permit  evaluation 
of  the  theoretical  statistical  moments  of  the  laser  beam  intensity  fluctua¬ 
tions,  its  functional  form  did  not  lend  itself  to  calculations  involving  the 
probability  of  detection  which  is  the  cumulative  fade-time.  So  as  to 
calculate  the  cumulative  fade  time  at  least  at  the  long  ranges,  the 
so-called  K-distribution  was  used.  The  data  from  the  March  1980 
experiments  did  show  that  the  K-distribution  was  a  reasonable  approxi¬ 
mation  to  the  laser  light  scattering  at  the  long  ranges.  It  was  the 
purpose  of  the  present  study  to  refine  the  earlier  proposed  universal 
math  model  so  that  the  probabilities  of  detection  for  various  thresholds 
could  be  calculated.  Additional  data  was  obtained  during  the  last  days 
in  August  1980  at  the  Space  Shuttle  Landing  Facility. 
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EXPERIMENTAL  MEASUREMENTS 

The  propagation  experiments  conducted  during  August  1980  were 
the  same  as  the  earlier  experiments  in  March  1980  and  they  are  described 
in  Reference  [1],  These  new  measurements  indicate  that  the  data  at  a 
particular  range  from  the  transmitter  can  vary  significantly  depending 
upon  the  strength  of  the  turbulence,  C*.  The  C*  depends  upon  the  air 
temperature,  time  of  day,  and  even  time  of  year.  As  in  the  previous 
data,  we  found  that  the  mean-square  fluctuations  of  the  laser  intensity 
normalized  with  respect  to  the  square  of  the  mean  value  of  the  laser 
signal  reached  a  maximum  value  of  approximately  6.  From  this  high 
value,  the  normalized  mean-square  fluctuations  decreased  to  the  limiting 
value  of  2. 

Due  to  the  strength  of  the  turbulence  during  the  August  measure¬ 
ments,  we  found  that  the  normalized  moments  of  the  laser  intensity 
fluctuations  at  ranges  of  2,500  -  3,000  meters  were  close  to  those  of  the 
negative  exponential  distribution.  The  negative  exponential  has  long 
been  recognized  theoretically  as  the  limiting  distribution  for  the  inten¬ 
sity  fluctuations  at  very  long  distances  under  conditions  of  very  strong 
turbulence.  We  believe  this  is  the  first  time  ever  any  research  team 
has  obtained  data  which  comes  so  close  to  supporting  this  particular 
mathematical  theory. 

RESULTS 

As  in  the  previous  report,  the  first  five  normalized  moments  <ln>/ 
<I>n,  n  =  1,2, 3, 4, 5  associated  with  the  measured  intensity  fluctuations 
of  the  optical  beam  were  calculated.  In  Figure  3-1  we  have  plotted  the 
third,  fourth,  and  fifth  normalized  moments  as  a  function  of  the  second 
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normalized  moment  <I2>/<I>2.  The  same  data  as  taken  in  March  is 

shown  in  Figure  3-2.  In  Figure  3-1  the  data  appears  to  be 

widely  dispersed  as  compared  with  Figure  3-2.  This  is  partly  due  to 

the  fact  that  there  was  a  lot  more  fluctuation  in  the  C2  from  minute  to 

n 

minute  during  the  August  measurements  as  compared  with  the  March 
measurements. 

To  better  understand  how  the  data  changes  as  a  function  of  dis¬ 
tance  and  strength  of  turbulence,  C2,we  have  plotted  the  third  normal¬ 
ized  moment  only  as  a  function  of  the  second  normalized  moment  in 
Figures  3-3,  and  3-4,  illustrating  the  data  at  various  distances  and  time 
of  the  day.  The  temperature  is  also  a  big  factor  in  this  data.  To  see 
how  the  data  moves  as  a  function  of  distance  and/or  structure  constant 
C2  ,  we  have  drawn  a  curve  that  we  believe  the  data  follows.  That  is 
to  say,  if  data  at  various  distances  could  all  be  obtained  simultaneously 
for  a  particular  value  of  the  structure  constant  C2,  we  believe  the 
data  would  move  along  the  lower  curve  shown  as  a  function  of  distance 
until  the  normalized  second  moment,  <I2>/<I>2,  reached  a  value  near  6 
and  then  as  the  distance  increased  even  more  the  data  would  be  dis¬ 
persed  along  the  upper  curve  approaching  the  limiting  point  where 
<I2>/<I>2  =  2.  Since  the  actual  data  obtained  was  gathered  at  different 
times  of  the  day,  the  value  of  the  structure  constant  C2  changed  by  as 
much  as  two  orders  of  magnitude  during  this  time  frame.  This  accounts 
for  the  fact  that  the  data  taken  in  August  at  500  meters,  which  was 
taken  from  2:15  p.m.  to  2:51  p.m.,  appears  to  be  closer  to  the  limiting 
value  <I2>/<I>2  =  2  than  does  the  data  at  1,000  meters  which  was  taken 
from  4:29  p.m.  to  4:56  p.m.  Similarly,  the  data  at  2,500  meters  was 
taken  from  2:32  p.m.  to  3:08  p.m.  while  the  data  at  3,000  meters  was 
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Fig.  3-4.  Measured  values  of  the  normalized  third  moment  obtained 
during  August,  1980  experiments.  The  curve  indicates  movement 
of  the  data  as  a  function  of  distance  and/or  structure  constant. 
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taken  from  4:13  p.m.  to  4:43  p.m.  Thus  we  were  actually  closer  to  the 
limiting  value  <I2>/<1>2  =  2  at  2,500  meters  than  we  were  at  3,000 
meters.  The  time  of  day  and  temperature  are  very  important  factors, 
therefore,  in  determining  the  strength  of  the  turbulence  as  well  as  the 
distance.  During  the  August  measurements,  the  value  of  the  structure 

2  .13  .2/3  .11  .2/3 

constant  C  varied  from  10  m  to  10  m  .  The  values 
n 

correspond  to  conditions  of  very  strong  turbulence. 

UNIVERSAL  MATH  MODEL 

Relying  on  arguments  similar  to  those  given  in  Section  V'  of  Refer¬ 
ence  (1],  we  assume  the  received  field  is  composed  of  two  principal  compo¬ 
nents  —  the  specular  component  resulting  from  the  forward  scattering 
by  the  large  eddies  along  the  propagation  axis  and  the  diffuse  component 
caused  by  multiple  scattering  by  the  off-axis  eddies.  Thus  the  received 
field  is  mathematically  represented  by 

U(t)  =  (Ae10  +  Relc()  eluit,  (3.1) 

where  the  first  term  corresponds  to  the  specular  component  and  the 
second  term  the  diffuse  component. 

Over  short  propagation  paths  the  field  U(t)  consists  primarily  of 
the  specular  component.  Many  arguments  have  been  presented  to 
suggest  that  the  amplitude  A  of  this  specular  component  satisfies  Log- 
Normal  statistics.  There  is  an  abundance  of  data  to  support  this  assump¬ 
tion,  including  our  own  short-range  data  obtained  during  March  and 
August  experiments  [1].  However,  over  longer  propagation  paths  it  is 
not  clear  what  statistics  the  amplitude  A  will  satisfy  since  the  amplitude 


25 


of  the  field  will  consist  of  both  specular  and  diffuse  components  rather 
than  the  specular  component  alone.  In  our  previous  model  we  made  the 
assumption  that  A  remained  Log-Normal  for  all  propagation  path  lengths. 
From  a  modeling  standpoint  this  assumption  limits  some  of  the  flexibility 
of  the  model  to  fit  the  experimental  data.  We  now  believe  that  it  is 
more  meaningful  to  assume  that  A  is  associated  with  a  more  general 
distribution  than  the  Log-Normal  distribution,  but  one  which  exhibits 
the  characteristics  of  the  Log-Normal  distribution  at  least  over  short 
propagation  path  lengths.  That  is  to  say,  the  distribution  should 
contain  certain  parameters  that  can  be  selected  from  comparisons  with 
experimental  data  and  which  have  the  effect  of  changing  the  character¬ 
istics  of  the  probability  distribution  as  the  parameters  change.  A 
general  distribution  of  this  type  is  the  m-distribution  defined  by 

p(A)  =  ,-MAVc,  A  >  0  (3.2) 

r(M)  cM 

where  c  =  <A2>  is  the  average  of  the  intensity  A2  and  M  is  the  primary 
parameter  that  controls  the  characteristics  of  the  distribution.  The 
function  T"(.)  is  the  gamma  function  ([3],  pp.  253-294). 

For  large  values  of  M,  the  m-distribution  behaves  like  a  Log- 
Normal  distribution  (see  Appendix  A),  for  M  =  \  it  is  Gaussian  in 
functional  form,  while  for  M  =  1  it  becomes  the  well-known  Rayleigh 
distribution.  Hence  the  m-distribution  seems  appropriate  here  because 
of  the  close  relationship  it  bears  to  the  standard  distributions  associated 
with  communication  theory. 
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The  amplitude  R  of  the  diffuse  component  is  also  assumed  to  be 
m-distributed,  i.e.,  / 


D2m'l  _D  2  /u 

p(R)  =  — — R  ■  e"  R  /b 


(3.3) 


F(m)  b 


m 


This  is  the  same  assumption  we  previously  made  in  Reference  [1]  and 
the  reasons  for  it  are  likewise  stated  there.  Also  as  before,  we  assume 
that  both  phase  angles  0  and  $  are  uniformly  distributed,  and  finally 
that  all  random  variables  are  independent. 


A.  Development  of  the  PDF : 

Our  previous  development  of  the  PDF  given  in  Reference  [1]  lead 
to  an  expression  that  was  difficult  to  use  for  numerical  calcalations  of 
the  cumulative  distribution  function.  The  present  development,  although 
similar  in  nature,  leads  to  a  simpler  functional  form  from  which  calcula¬ 
tions  can  be  more  easily  obtained.  The  intensity  of  the  field  (3.1)  is 

I  =  |U(t)  |2  =  A2  +  R2  +  2AR  cos  (0-0).  (3.4) 


Following  the  technique  outlined  in  [1],  the  PDF  is  determined  from  the 
relation 


P(I)  = 


where 


a 


z  C(z)  J0  (Vi  z)  dz 


/ 


0 

oo 


C(z)  =  /  p(A)  JAAz)  dA  /  p(R)  J_.(Rz)  dR. 


00 

/p 


(3.5) 


(3.6) 


The  function  J0  (.)  is  the  standard  Bessel  function  of  zero  order  ([3], 
pp.  355-434).  The  evaluation  of  these  integrals,  which  is  given  in 
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Appendix  B,  leads  to  the  infinite  series  representation 


P(I)  = 


mr(m) 

EfXM) 


e-ml/b 


k=0 


(3.7) 

r(M  +  i)  (rmW 
r (m  -  k  +  j) 


where 

k\  =  k! 

i/  ■  i)!  (3.8) 

is  the  binomial  coefficient  and  denotes  the  Laguerre  polynomials  [3]. 

The  p aramter  r  appearing  in  (3.7)  is  the  power  ratio  of  the  spec¬ 
ular  to  diffuse  components,  i.e., 

r  =  c/b.  (3.9) 

Over  short  propagation  path  lengths  the  parameter  r  is  very  large  and 
monotomically  decreases  to  zero  as  the  strength  of  the  turbulence  in¬ 
creases  and  the  path  length  increases.  In  fact,  the  limiting  form  of 
(3.7)  is  obtained  by  letting  r  -»  0  and  m  1 ,  which  leads  to 


PCD  ■  i  e-‘/b 


<3 


This  limiting  form  is  the  negative  exponential  distribution  which  has 
been  strongly  supported  by  reasearchers  over  the  years. 

B .  Theoretical  Moments: 

The  statistical  moments  of  the  intensity  1  are  computed  from  the 
expression 

co 

<In>  =  j  In  p(I)  dl,  n  =  1,2,...  (3.11) 


.10) 
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from  which  we  get  (see  Appendix  C) 


n 


k=Q 


r(M  +  k)  r(m  +  n  -  k)  /rm\k 
r(M)  r (m)  \M  ) 


(3.12) 


For  purposes  of  comparing  the  theoretical  moments  with  those  obtained 
from  experimental  data  it  is  necessary  to  normalize  them.  This  we  do 
by  first  setting  n=l  to  find  the  mean 


<I>  =  b(l  +  r), 


(3.13) 


and  then  dividing  (3.12)  by  <I>n. 

n 


This  action  results  in  the  expression 


<In>  =  1  /n\2 

<I>n  (1  +  r)n  W 

k=0 

where  we  have  introduced  the  notation 


Mn-k  ak 


(3.14) 


..  _  <RZ1S  _  F(m  +  k) 
Mk  -  IT 

K  <R^>k  mK  T(m) 

and 

_ <A2k>  _  r(M  +  k) 

ak~  1  Tc  -  H - 

<A  >  Mk  r(M) 


(3.15) 

(3.16) 


Eq.  (3.14)  has  the  same  functional  form  as  Eq.  (5.19)  in  Reference 
[1].  The  primary  distinction  is  that  the  normalized  moments  ak  of  the 
specular  component  and  now  associated  with  the  m-distribution  rather 
than  the  Log-Normal  distribution. 

In  order  to  match  the  theoretical  moments  (3.14)  with  actual  data 
obtained  from  field  experiments,  we  have  to  appropriately  select  three 
parameters:  m,  M  and  r.  We  do  this  in  such  a  way  that  the  first  two 
normalized  moments  of  I  match  the  data  identically  and  the  third  moment 
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Fig.  3-5.  Parameters  of  the  m-distributions  associated  vith  the 

diffuse  and  specular  components  of  the  laser  beam  as  a  function 
of  the  normalized  variance. 
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3-6.  Measured  values  of  the  normalized  third,  fourth  and  fifth 
moments  obtained  during  March,  1980  experiments  and  compared  wit) 
values  expected  from  the  universal  model. 


Distance  Time  of  Day 
-  500  m  (2:15  p.m.-2:51 
”1000  m  (4:29  p.m.-4:56 


urvc 


is  then  matched  by  some  "best  fit"  criterion. 

The  parameters  m  and  M  as  a  function  of  the  normalized  variance 
Oj2  that  appeared  to  give  the  best  visual  fit  to  the  data  are  graphed  in 
Figure  3-5.  The  remaining  parameter  r  is  then  calculated  from  the 
formula 


r  =  2  -  x  +  V  (a2  +  M2  EjQ  X  I  4  -,a^2  (3>17 

X  ~  ^2 

where  a2  =  1  +  1/M,  \i2  =  1  +  1/m  and  x  =  <I2>/<I>2.  Eq.  (3.17)  is 
equivalent  to  Eq.  (5.24)  in  [1].  The  resulting  theoretical  curves  along 
with  the  data  obtained  in  March  and  August,  respectively,  are  shown  in 
Figures  3-6  and  3-7. 

The  three  parameters  m,  M  and  r  of  our  theoretical  model  have 
been  tied  to  the  normalized  variance  Oj2  as  suggested  above.  The  vari¬ 
ance  as  a  function  of  the  parameter  =  (0.5  C2  k7^6  L11^6)'5,  where 
C2  is  the  structure  constant,  k  is  the  wave  number  of  the  optical 
beam,  and  L  is  the  path  length,  is  shown  in  Figure  3-8  for  the  actual 
data  obtained  in  August.  The  data  is  widely  scattered,  as  discovered 
also  by  other  researchers,  so  that  a  single  curve  drawn  through  the 
data  does  not  seem  feasible  at  this  time.  The  dashed  curves  that  we 
have  drawn  suggest  possible  upper  and  lower  bounds  for  the  normalized 
variance  as  a  function  of  PQ. 

C.  Cumulative  Distribution: 

To  derive  the  cumulative  distribution  function  from  which  the 
fading  probabilities  (and  hence,  the  detection  probabilities)  can  be 
computed,  we  first  normalize  the  density  function  (3.7)  so  that  the 
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mean  value  is  unity.  Thus  by  introducing  the  normalized  variable  y  = 
!/<!>,  we  see  that  (3.7)  becomes 


p(y)  =  (1  *  r>  mr(m?  e-«  kL.  1(1  ♦  r)  my] 

r(M)  ^  k!  K 

k=0 

k 

X  "WiT  r(M  v  i)  (rm/M)j 

Z-A]/  r(m  +  j  -  k) 

j=o 

The  probability  of  detection  is  then  defined  by 

pd  =  J  p(y)  dy 

T 

for  a  specified  threshold  value  T  =  IT/I0,  where  I  is  the  mean 
intensity.  Through  application  of  the  integral  formula  ([4],  p. 


00 

J e"x  Ln(x)  dx  =  e"y  (Ln(y)  -  Ln_j(y)] 

y 

we  immediately  find  that 


£ 


(-1^  [Lk(z)  -  L.  .(z)] 

TT  K  K  1 


X 


j=o 


r(M  ♦  i) 
r(m  +T-  k) 


(rm/M)^ 


) 


where 


z  =  (1  +  r)  mT. 


The  probability  of  fade  is  therefore 


P 


f 


=  1 


P 


d 


(3.18) 

(3.19) 

value  of 
844) 

(3.20) 


(3.21 ) 

(3.22) 

(3.23) 
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The  fade  and  detection  probabilities  associated  with  long-range 

propagation  paths  for  which  Oj2  =  1.5  are  shown  in  Figure  3-9  for  the 

universal  math  model  described  here  as  well  as  similar  curves  predicted 

by  the  Log-Normal  and  K-distribution  models.  For  low  threshold  values 

we  find  the  Log-Normal  model  predicting  a  much  higher  probability  of 

detection  than  our  universal  model,  while  the  K-distribution  gives 

cumulative  probabilities  that  closely  resemble  the  universal  math  model 

as  we  previously  suggested  in  Reference  [1],  As  the  variance  de- 

2 

creases  for  longer  ranges  to  the  limit  value  Oj  =  1 ,  we  have  found  that 

the  K-distribution  curves  approach  those  predicted  by  the  negative 

exponential  model  in  a  manner  similar  to  our  universal  model.  In  Tigure 

3-10  we  have  shown  how  this  limiting  cumulative  distribution  curve 

(depicted  by  the  dashed  curve)  is  approached  for  decreasing  variance 

2 

(corresponding  to  increased  path  lengths)  from  a  value  of  Oj  =  3.5  to 


2 


2 

For  variances  such  that  Oj  >  1.5  at  long  distances  the  K-distribution 
deviates  more  strongly  from  the  cumulative  distribution  curves  predicted 
by  the  universal  math  model.  For  that  reason  we  do  not  believe  the 
K-distribution  represents  a  good  model  for  these  values  of  the  variance 
which  correspond  to  more  intermediate  path  lengths  and  conditions  of 
turbulence. 

For  short  ranges  and  thresholds  such  that  T  £  0.4  we  find  little 
distinction  between  the  cumulative  probability  curves  as  predicted  by 
the  Log-Normal  model  and  the  universal  model.  However,  for  lower 
threshold  values  the  log-normal  model  predicts  a  higher  probability  of 
detection  than  does  the  universal  model.  Some  representative  curves  of 
the  universal  model  for  the  cumulative  probabilities  are  shown  in  Figure 
3-11  corresponding  to  variances  of  Oj2  =  0.5,  1.0,  2.5. 
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PROB.  OF  FADE 


THRESHOLD  (  /T  / 1()) 

Fig.  3-9.  Comparison  of  cumulative  fading  models  at  long  ranges 
for  the  K-distribution ,  the  Log-Normal  distribution  and  the 
universal  model  as  a  function  of  the  ratio  of  receiver  threshold 
to  mean  intensity. 
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PROB  OF  DETECTION 


PROB.  OF  FADE 


PROB  OF  DETECTION 


Fig.  3-11.  Probability  of  fade  curves  as  predicted  by  the  universal 
model  for  short  ranges  with  different  variances. 
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SECTION  IV 


CONCLUSIONS 


The  MILES  system  was  designed  as  a  tactical  weapon  fire  simulator. 

The  kill  zone  was  tailored  to  have  a  shape  so  as  to  have  an  approximate 
kill  probability  equal  to  the  kill  probability  of  a  particular  weapon 
engaging  a  particular  target.  To  do  this  the  beam  was  made  wide 
enough  to  illuminate  several  detectors  for  any  well-aimed  shot.  In  some 
instances,  the  entire  target  would  be  illuminated  by  the  laser  beam. 

The  output  of  each  of  the  illuminated  detectors  was  summed  and  then  discrim¬ 
inated  against  a  fixed  threshold.  There  was  no  attempt  to  simulate  the 
hit  of  a  round  on  any  particular  location  on  the  target.  The  tested 
MILES  system  as  currently  configured  would  make  a  very  poor  gunnery 
marksmanship  trainer.  The  following  list  of  recommendations  are  pro¬ 
posed  as  MILES  add-on  which  would  allow  the  MILES  system  to  be  used 
as  a  marksmanship  trainer. 

RECOMMENDATIONS 

1.  Replace  MILES  solar  cell  detectors  with  reversed  biased  photo¬ 
diodes. 

2.  Replace  MILES  detector  belts  with  a  matrix  of  independent  detec¬ 
tors  each  with  its  own  pre-amplifier  and  discriminator  threshold. 

The  MILES  pre-amplifier  and  discriminator  can  be  retained. 

3.  Retain  the  MILES  decoder  box,  but  one  must  be  used  for  each 
detector. 
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4.  Replace  MILES  laser  transmitter  with  transmitter  with  much  smaller 
beam  divergence  and  more  beam  uniformity.  This  will  ensure  a 
single  narrow  gaussian  beam  as  opposed  to  the  wide  MILES  tri-gaussian 
beam. 

5.  Retain  the  MILES  encoder  which  is  used  to  drive  the  transmitter 
with  standard  MILES  codes. 

If  these  recommendations  are  implemented,  the  following  questions  would 
require  answers: 

a.  What  should  be  the  transmitter  power? 

b.  What  should  be  the  transmitter  beam  divergence? 

c.  What  should  be  the  photo-diode  detector  spacing  on  the 

target? 

To  begin  to  answer  these  system  design  questions,  the  receiver  struc¬ 
ture  must  be  defined  and  the  atmospheric  optical  communication  channel 
statistics  must  be  modelled.  The  NTEC  research  team  proposed  a  CW 
laser  transmitter  operating  at  a  high  PRF  and  a  phase-locked  loop  or  a 
matched  filter  receiver.  The  NTEC  receivers  and  CW  laser  high  PRF 
transmitter  also  requires  a  knowledge  of  the  optical  communication 
channel  statistics  so  that  the  pulse  rate  of  the  transmitter  can  be  set 
and  the  optimum  threshold  of  the  receiver  can  be  established. 

The  statistical  model  developed  by  the  UCF  research  team  estab¬ 
lished  the  statistical  moments  and  the  cummulative  fading  as  a  function 
of  the  threshold  setting  of  the  discriminator.  The  following  observa¬ 
tions  can  be  made  about  the  statistical  model  of  the  optical  communica¬ 
tion  channel: 
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1.  The  measured  long  range  statistics  were  Negative  Exponential  as 
shown  in  Equation  3-10.  This  is  the  first  time  this  statistic  has 
ever  been  measured  although  it  has  been  theoretically  predicted 
for  some  ten  years. 

2.  The  developed  math  model  matches  the  Log-Normal  statistic  at 
short  ranges  and  the  asymptotic  long-range  Negative  Exponential 
statistic . 

3.  The  intermediate  ranges  have  the  lowest  probability  of  detection  or 
rather  the  greatest  cummulative  fading. 

4.  The  detection  probability  gets  better  at  the  very  long  ranges. 

5.  The  amount  of  scattering  can  vary  considerably  at  a  fixed  range. 
For  example  at  2,000  meters  the  detection  probability  may  be  at 
one  value  and  then  when  the  scattering  gets  stronger  the  detection 
probability  may  actually  improve.  Conversely,  if  the  scattering 
would  weaken,  it  is  possible  the  detection  probability  may  actually 
get  worse. 

6.  The  threshold  setting  of  the  receiver  should  be  determined  using 
tne  worst  case  using  intermediate  range  statistics. 

The  system  parameters  can  be  set  using  the  math  model  developed 
here,  but  to  set  the  PRF  of  the  proposed  NTEC  CW  laser  transmitter, 
the  average  length  of  fade  time  needs  to  be  determined.  This  math 
model  was  not  part  of  this  project,  but  if  the  development  of  the  pro¬ 
posed  NTEC  system  continues,  that  math  model  must  be  constructed. 

The  m-distribution  used  in  the  analysis  would  form  the  basis  for  the 
math  model  of  the  average  fade  period  as  a  function  of  the  receiver 
threshold. 
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The  combined  UCF  and  NTEC  study  suggests  that  the  CV\  laser 
high  PRF  transmitter  and  receiver  be  developed  as  MILES  add-ons  for 
long-range  gunnery  marksmanship  training.  To  set  the  parameters  of 
that  CW  laser  system,  a  math  model  for  the  average  fade  period  would 
be  required. 
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APPENDIX  A 


LOG-NORMAL  APPROXIMATION  TO  THE  m-DISTRIBUTlON 


The  m-distribution 


P(A)  = 


2M 


M 


T(M)c 


M 


a2M-1  e-MA2/c 


( A  - 1  ) 


for  large  values  of  the  parameter  M,  can  be  closely  approximated  by  the 
Log-Normal  distribution.  To  see  how  this  comes  about,  we  start  by  intro¬ 
ducing  the  change  of  variables 


A2  _ 


c 


(A-2 ) 


so  that  (A-l)  becomes 


P(x)  = 


M 

7TM) 


eMx  exp(-Mex). 


(A-3) 


Under  the  assumption  that  x  <  <  1  (i.e.,  A2  <  <  c),  we  can  use  the 

approximation 

ex  =  1  +  x  +  \  x2 

in  (A-3)  from  which  we  get 


M 

p(x)  -  TIM)  ^  expl'M(1  +  x  +  15x2)1 


iW 


~  — —  e'M  exp(-Mx2/2). 


( A-4) 


Now  returning  to  the  original  variable  by  use  of  (A-2),  we  see  that 


omM  -M 

P<A>  s  TTMTA"  exp 


^  ^log  A2  -  log  c^z 


(A-5) 


Finally,  we  invoke  Stirling's  formula  (13],  p.  257) 


-M 


r(M)  -  e  M  ,  M  >  >  1 


from  which  we  deduce 


p(A)  s  exp 


Pm/  2  \2 

y  log  A  -  log  c 


or  equivalently 


P(A)  S 


v2rt  A 


exp 


-(log  A  -  h  log  c) 
2  °LN2 


(A-6) 


where  we  have  set  M  =  1  *  ^earty'  the  Log-Normal 


distribution. 


APPENDIX  B 


DERIVATION  OF  THE  PDF 


From  Eq.  (5.6)  in  Reference  [1],  we  have 


/ 


00 


p(R)  J0(Rz)  dR  =  1F1(m;  1;  -bz2/4m) 


(  B  - 1 ) 


where  iF^a;  b;  c)  is  the  confluent  hypergeometric  function  ([3],  pp. 


503-536).  Similarly  we  find 


/ 


00 


p(A)  J0(Az)  dz  =  iFj(M;  1;  -  cz2/4M) 


so  that  Eq.  (3.6)  in  Section  III  leads  to  the  result 


(B-2 ) 


C(z)  =  1F1(M;  1;  -  cz2/4M)  jFjfm;  1;  -bz2/4m) 


(B-3) 


Next  we  employ  Rummer's  transformation  ([3],  p.  505) 


iF j(a;  b;  -x)  =  e’x  jF^b  -a;  b;  x) 


so  that 


C(z) 


=  e"bz2/4m  iF,(M;  1;  -  cz2/4M)  1F1(1  -  m;  1;  bz2/4m),  (B-4) 


and  finally  we  use  the  formula  ((5],  p.  187) 


jFjfa;  c;  px)  iFi(a’;  c';  qx) 


00 


Z(a) k (px)  3F2(a',  1  -  c  -  k,  -  k;  c',  1  -  a  -  k;  -  q/p) 
k!  (c). 


k=0 


where 


-  r(a  +  k  )  ,  k  =  0,1,2, 
c  '  k"  r(a) 
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is  the  Pochhammer  symbol  and  3F2(a,b,c;  d,e;  x)  is  a  generalized 
hypergeometric  function,  to  obtain 


00 


C(z)  =  e 


-bz2/4m 


r 


k=0 


(1  -  m)k  (bz2/4m)k 
(k!)Z 


(B-5) 


x  3P 2  (M ,  -k,  -k;  1,  m  -  k;  cm/bM). 


The  substitution  of  (B-5)  into  the  integral 


P(I) 


00 

*5  J  ZC 


zC(z)  Jq(JTz)  dz 


(B-6) 


and  the  subsequent  interchange  of  integration  with  summation,  together 
with  the  integral  formula  ([4],  p.  718) 


/' 


-x2  ..2n+l 


e"A  x"“Ti  Jq(2x  Vz)  dx  =  e”z  Ln(z), 


leads  to  the  expression 


00 


PCD  =  g  emI/b 


Z(1  -  m)k  Lk(ml/b) 
R1 - 

k=0 


(B-7) 


x  3B2(M,  -k,  -k;  1,  m  -  k;  cm/bM) 


Replacing  3F2(;  ;)  with  its  series  representation  and  simplifying,  we 
finally  get  Eq.  (3-7). 
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APPENDIX  C 


DERIVATION  OF  THE  THEORETICAL  MOMENTS 


Starting  with  the  equation 


in> .  / 


P(U  dl,  n  =  1,2 . 


(C-l) 


where  p(I)  is  given  by  (3.7),  and  employing  the  integral  formula  ([4], 


p.  844) 


/• 

e'st  tP  Ln(t)  dt  =  r(p  +  1)  s"(P  +  !)  F(-n,  p  +  1;  1;  1/s) 


we  find 


<rn>  -  r<m)  V*  (-D  n!  F(’k'  n  +  1;  1;  1) 

\m)  rm  irr 


E  ©'  „Kt4 


)  (rm/M)^, 


(C-2) 


where  F(a,b;  c;  x)  is  the  hypergeometric  function.  Using  the  relation 
([31,  P-  556) 


F(a,b;  c;  1)  =  F(c)  T(c  -  a  -  b) 
F(c  -  a)  F(c  -  b) 


we  find  in  our  case 


F(-k,  n  +  1;  1;  1)  = 


.  l(-l)k  n!/k!(n  -  k)f  ,  kin 


,  k  >  n  (C-3) 


and  therefore  (C-2)  becomes 


(C-4) 


<>">  4 


n 


_ (n!)2  l~(m)  r(M  +  j)  (rm/M)j 

(n  -  k) !  (j !  )z  [  (k  -  ITTp  r(M)  T(m  +  j  -  k) 


k=0  j=0 


Eq.  (C-4)  has  the  functional  form 


<>">  ■  (a 


n 


n  k 


C(k,  j) 


k=0  j-0 

so  that  making  the  chang j  of  indices  p  =  k  -  j  we  have 


n  n-j 


<'">  *  a 


£  23 c<p  *  >■ » 

j=0  p=o 


or 


(C-5 ) 


rn,  _  fb\n  ^  (n!)2  I~(M  +  j)  (rm/M)] 


<IU>  = 


m /  £_^  r(M)  (j!)2  (n  -  j)! 

j=0 
n-j 


v  V  (n  -  j)!  r(m) 

x  /  v  (n  -  p  -  j)!  r(m  -  p)(p!)2 

p=0 

However,  the  second  series  in  (C-6)  reduces  to 


n-J 

£ 

p=0 


_ (n  -  j)!  r(m) 

(n  -  p  -  j)!  r(m  -  p)  (pTF 


n-j 

£ 

p=0 


(1  -  m)  (-n  *  J)  0J> 
(T7p -  pT 


=  F(1  -  m,  -n  +  j;  1;  1) 


f(m  +  n  -  j)  ,  (C-7 ) 

r(m)  T(n  -  j  +  1) 
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and  substituting  this  expression  into  (C-6)  gives  the  final  result 
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Modeling  solar  cells  for  use  as  optical  detectors: 
background  illumination  effects 

Leo  A.  Mallette  and  Ronald  L.  Phillips 

Solar  cells  have  traditionally  been  used  lor  direct  sunlight  to  energy  conversion,  hut  there  ha-  been  rt  l.it i\  eh 
little  investigation  into  their  use  as  a  low  data  rate  optical  detector  This  paper  presents  an  experiment. i! 
procedure  used  to  determine  the  ac  model  of  a  specific  solar  cell  A  lumped  circuit  model  and  governing 
equations  are  developed.  Open  circuit  responses  to  pulses  are  used  to  determine  values  lor  the  internal 
capacitances  as  a  function  of  background  illumination  Possible  problems  encountered  with  noise  genera 
tion  are  also  reviewed. 


Introduction 

List  of  Symbols 

/  terminal  current  into  solar  cell. 
h  light  generated  current  internal  to  cell, 
kd  constant  coefficient  associated  with  fa, 

k,  constant  coefficient  associated  with 

V  terminal  voltage  of  cell. 

I'd,  diffusion  or  built  in  potential. 

rd  time  constant  associated  with  both  Cd  and 
r,  time  constant  associated  with  and 
tp  lifetime  of  minority  carriers. 

Genera! 

The  photovoltaic  effect  was  first  observed  by  Bec- 
querel  in  1839,  and  over  100  yr  later,  it  was  observed  by 
Ohl  in  a  silicon  p-n  junction.1  Modern  manufacturing 
techniques  have  made  low-cost  silicon  solar  cells  readily 
available.  Availability  of  the  material  and  price  have 
prompted  investigation  into  their  use  in  areas  other 
than  direct  sunlight  to  energy  conversion;  however,  the 
feasibility  of  using  a  solar  cell  as  a  low  data  rate  optical 
detector  has  been  investigated  by  relatively  few  re¬ 
searchers.2-1  In  order  to  predict  the  response  of  a  solar 
cell,  the  transfer  function  or  the  equivalent  circuit  must 
be  known.  In  this  paper  an  earlier  model1  is  assumed, 
the  individual  elements  are  explained,  the  governing 
equations  are  given,  and  a  method  to  fit  numbers  to  the 
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equations  is  developed.  This  paper  assumes  a  constant 
temperature  of  300  K;  variation  with  temperature  will 
be  presented  in  a  later  paper. 

Model 

A  solar  cell  is  a  p-n  junction  with  part  of  it  surf  ace 
exposed  to  the  environment.  A  diode  and  current 
generator  in  parallel  can  be  used  as  a  rough  approxi¬ 
mation  for  a  solar  cell  model.  All  diodes  have  an  in¬ 
herent  series  resistance  Rs  and  shunt  resistance  Rsn-: 
and  incorporating  these  yields  an  excellent  approxi¬ 
mation  to  the  dc  model  of  the  solar  cell.  The  cell  cur¬ 
rent  is  proportional  to  the  incident  power  by  the  sped  ral 
responsivity  constant  and  is  given  by4 

/  =  -/»  +  / n  |exp|-q(  V  -  IRsl'ik  7j  -  11 

+  ( V  -  IRsI/Rsh  lamperesi  ill 

For  practical  purposes,  Shirland4  defined  Rs  and  Rsh 
in  terms  of  the  solar  cell  volt-ampere  ( \*-/l  character¬ 
istics  as  follows: 

Rs  is  the  reciprocal  of  the  slope  of  the 

Y-I  curve  at  /  =  0  (amperes',  rji 

Rsh  is  the  reciprocal  of  the  slope  of  the 

V  -/ curve  at  V  *=  o  (volto  <;<» 

The  remaining  unknown  element  in  the  model  (Fig.  1 ) 
is  the  diode,  which  can  be  characterized  by  three  term-': 
(a)  the  junction  resistance  r,  ( b )  the  junction  capaci¬ 
tance  Cj,  and  (c)  the  diffusion  capacitance  C,/. 

The  junction  resistance,  usually  called  the  dynamic 
resistance,  is  defined  as  the  reciprocal  of  the  slope  of  the 
V-I  curve  and  is  given  bvr’ 

rikT 

r  • - exp(-<)V/)/i|*7')  (ohms)  < 4 ) 

<?/ n 
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Fig  2.  Volt-ampere  characteristics  for  photovoltaic  cell  with  and 
without  background  illumination. 


The  junction  capacitance  is  defined'5  as  the  incremental 
charge  increase  dq  of  fixed,  unneutralized  donors  or 
acceptors  per  incremental  voltage  change  dt  and  is 
given  by’ 7  81 


Cj  -  V<V»,  -  '7)>'  '  (faradsl.  lot 

The  diffusion  capacitance  term  C'd  is  a  result  of  the 
photoinjected  electron-hole  pairs.  On  either  side  of  the 
junction,  pairs  of  majority  and  minority  carriers  are 
generated  whenever  a  photon  is  absorbed.  The  addi¬ 
tion  of  majority  carriers  is  negligible,  but  there  is  a 
highly  noticeable  increase  in  the  percentage  of  minority 
carriers.  These  carriers  will  either  recombine  with 
majority  carriers  or  diffuse  to  the  junction  and  drift  to 
the  opposite  side  of  the  junction  where  they  become 
majority  carriers.  The  lifetime  of  the  minority  carrier 
is  defined  as  rp,6  and  the  diffusion  capacitance  can  be 
shown  to  be* 


Cd  =  rp/r  (farads).  (fi) 

where  r  is  defined  by  Eq.  (4). 

Since  rp  is  a  constant,  Cd  can  be  shown  to  be8 

C d  ~  kd  txpioV i>)  (farads)  (7) 

Experimental  Procedure 

Resistances 

Two  V-I  characteristics  of  the  solar  cell,  shown  in  Fig. 
2,  are  measured  using  a  commercially  available  curve 


tracer.  The  reciprocal  of  the  slope  of  the  V  I  curve, 
measured  in  the  presence  of  background  light  and 
evaluated  at  points  .4  (/-axis  intercept  i  and  H  i  I'-axi- 
intercept),  yields  the  shunt  and  series  resistances,  re 
spectively.4  This  curve  is  similar  to  the  no-light  curve, 
but  translated  down  by  /A  —  ( V  —  IRS)/R>H  amperes, 
and  translated  to  the  right  by  1RS  volts.  Note-  these 
V-I  characteristics  are  for  the  entire  solar  cell  and  not 
for  the  internal  diode. 

The  internal  diode  current  and  voltage  are  given 

by1 

V[,  =  \'  —  IHs.  (vohto.  i  **  i 

and 

Ip  =  I  —  i  V  —  //fs  )  ffs//  (amperes)  is 

The  data  points  of  the  V-I  curve  can  be  corrected  using 
Eqs.  (8)  and  (9)  to  supply  the  \'/,  -  In  curve  ol  the 
diode.  These  corrected  data  points  are  fitted  to  the 
equation* 

Id  =  /utexpiti \'h )  —  1 1  (ampere^  ■]". 

The  dynamic  resistance  r  is  obtained  by  taking  tin- 
reciprocal  of  the  derivative  of  Eq.  ( 101.  or  via  Eq.  1 4  i. 
noting  that  o  of  Eq.  (10)  is  q/qkT. 

Junction  Capacitance 

Saltsman*  used  the  circuit  of  Fig.  8  to  determine  the 
junction  capacitance  by  measuring  the  impulse  response 
of  the  cell  for  various  values  of  dc  bias.  The  bias  resistor 
Rh  was  chosen  so  that  Rs  «  Rh  «  Rsh  R,  and  ('. 
were  chosen  to  provide  ac  coupling,  to  keep  the  peak 
value  of  the  impulse  small,  and  to  have  a  long  time 
constant.  The  cell  voltage  V.  as  observed  on  an  oscil 
loscope.  is  a  decaying  exponential  superposed  on  the  dc 
bias,  which  is  generated  by  the  background  illumina¬ 
tion.  The  decay  time  constant  r,  will  vary  as  a  function 
of  temperature  and  background  illumination.  Tin- 
decay  time  constant  is  given  by* 

7.  J  HfA'  (set  I  till 

The  decay  time  constant  can  be  determined  from  the 
cell  voltage  V  as  observed  on  an  oscilloscope. 


cc 


Fig  .1  Test  circuit  used  to  determine  junction  cap.tt  il.tm «  «*t  tin 
solar  cell 
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Dividing  r;  by  Rp  yields  the  junction  capacitance  for 
various  values  of  reverse  bias  voltage.  These  values  of 
capacitance  must  be  fitted  to  Eq.  (5).  One  method  is 
to  assume  a  value  for  ->  and  plot  the  experimental  data 
points  as  shown  in  Fig.  4.  If  the  value  for  i  was  cor¬ 
rectly  chosen,  the  points  will  form  a  straight  line;  oth¬ 
erwise,  a  different  value  of  ■>  must  be  chosen.9  The  line 
connecting  the  experimentally  measured  points  may  be 
extended  to  the  V'/>  axis.  The  built-in  potential  can  be 
read  off  the  graph,  since  the  l’o-axis  intercept  is  at  ( Vb, 
—  \’[,\  =  0.  Similarly,  the  constant  coefficient  kj  can 
be  determined  from  the  l/O-axis  intercept. 

Diffusion  Capacitance 

The  authors  developed  a  method8  to  determine  the 
diffusion  capacitance  using  a  narrow-pulse  laser  and  a 
source  of  controllable  background  illumination.  The 
background  illumination  generates  a  dc  current  /*  and 
hence  a  voltage  Vp  across  the  diode.  The  two  lasers111 
coimpinged  on  the  face  of  the  solar  cell,  and  the  open 
circuit  decay  time  constant  was  measured  as  well  as  the 
dc  output  voltage.  The  magnitude  of  the  pulse  had  to 
be  kept  small,  compared  with  the  background  illumi¬ 
nation,  in  order  that  it  not  significantly  perturb  the  cell 
by  itself.  Due  to  background  illumination,  the  time 
constant  will  be  a  function  of  both  the  junction  capac¬ 
itance  and  the  diffusion  capacitance,  and  it  is  given 
by 

|  -  ‘  -  +  — +  -1  (Cj  +  Cd)  (seel.  (121 

L<  fis  +  Rj.)  Rsh  rJ 

The  l/(Rs  +  Rp)  term  is  negligible  because  of  the 
method  of  open  circuit  measurement  (coupling  directly 
to  an  oscilloscope).  The  other  terms  in  Eq.  (12)  have 
been  previously  determined  and  are  either  constants  or 
functions  of  Vp.11  A  graph  of  rj  vs  output  voltage  will 
yield  a  family  of  monotonically  decreasing  functions 
from  V  *  0  volts  to  V  *  V*,.  The  junction  capacitance 
dominates  at  the  lower  voltages,  and  (as  shown  below) 
noise  dominates  at  the  higher  voltages.  Since  kj  is  the 
only  unknown  in  Eq.  (12),  one  point  will  uniquely  de¬ 
termine  its  value.  Theoretically  this  is  sufficient,  but 
the  authors  have  found  that  fitting  an  entire  set  of 
measurements  will  eliminate  random  errors. 


Noise  Generation 

Because  of  internal  resistance,  physical  size,  exposure 
to  the  environment,  and  the  quantum  mechanism  <>t 
operation,  solar  cells  are  subject  to  nearly  every  known 
noise  identified  to  date!  In  considering  the  S.N'H.  -ix 
noises  the  reader  should  be  aware  of  are8-1* 

(a)  Johnson  noise.  Also  called  Nyquist.  resistance, 
or  thermal  noise. 

(b)  Temperature  noise  shifts  the  spectral  response 
toward  the  ir  as  temperature  increases. 

(c)  Modulation  noise.  Also  called  excess  or  i  / 
noise. 

(d)  Generation-recombination  (G-R)  noise.  The 
variation  can  be  caused  by  the  random  arrival  of  pho¬ 
tons  from  the  background  illumination.  For  this  case, 
the  noise  is  also  called  photon,  radiation,  or  background 
noise. 

(e)  Shot  noise  is  due  to  the  discrete  nature  of  the 
current.  The  total  current  of  the  unbiased  cell  is 
comprised  of  minute  current  pulses  produced  by  tilt- 
individual  electron-hole  pairs  being  generated.  Thi» 
noise  was  intense  when  the  cell  was  exposed  to  high 
levels  of  background  illumination. 

(f)  60-Hz  noise  is  the  modulation  of  the  background 
illumination  by  60  Hz  and  its  harmonics. 

Summary 

A  general  model  for  a  solar  cell  was  developed  using 
lumped  resistors  and  capacitors.  The  resistances  art- 
determined  from  the  V-I  characteristics  of  the  cell 
The  junction  and  diffusion  capacitances  are  determined 
by  observing  the  decay  time  constant  of  electrical  and 
optical  impulses,  respectively.  The  method  should  he 
applied  in  the  same  chronological  order  as  this  paper, 
isolating  the  ac  and  optical  effects,  then  releasing  the 
ac  isolation  requirement  to  measure  the  junction  ca¬ 
pacitance,  and  finally  the  optical  isolation  requirement 
to  calculate  the  diffusion  capacitance.  Changing  the 
background  illumination  will  cause  significant  variation 
in  the  value  of  r  and  hence  may  wreak  havoc  on  any 
form  of  digital  signal  processing  equipment  in  the  de¬ 
tection  system. 

The  authors  thank  W.  A.  Kownacki  for  his  many 
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